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Abstract We discuss the dynamics of a charged non- 
relativistic particle in electromagnetic field of a rotating 
magnetized celestial body. The equations of motion of 
the particle are obtained and some particular solutions 
are found. Effective potential energy is defined on the 
base of the first constant of motion. Regions accessi- 
ble and inaccessible for a charged particle motion are 
studied and depicted for different values of a constant 
of motion. 
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1 Introduction 

The field of magnetic dipole and motion of charged 
particles in this field have large practical significance 
for astrophysics. Magnetic field of many planets and 
stars can be thought of as a dipole field in good ap- 
proximation. Stationary case, when the magnetic mo- 
ment of a celestial body coincides with the rotation 
axis, is well studied. Particularly, motion of a charged 
particle in the Earth magnetic field is studied in de- 
tail. Solution of equation of motion for charged particle 
in a dipolar magnetic field gives trapping regions for 
particles, which have energy of limited range. These 



regions of planets are named radiation belts (Alfven 



1950 Holmes-Siedle & Adams 2002). The first theo- 



retical studies on the properties of the trajectories of 
a charged particle in the dipolar field where done by 



St0rmer (1907), DeVogelaere (1958) and Dragt (1965) 
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Rather complicated mathematical methods were devel- 
oped in order to calculate the trajectory of charged par- 



zano 



ticle in the dipole magnetic field (see for example ( Lan- 
|1968 )). Solution of the St0rmer's problem was ex- 



tended to magnetic field that is a superposition of the 
field of a dipole and a uniform magnetic field. Allowed 
and forbidden regions of the motion of charged parti- 
cles in such field was studied by |Katsiaris fe Psillakis| 
( fl~986 >. 

There are also well known celestial bodies, which di- 
rection of magnetic moment differs from the direction of 
axis of rotation. In this case electric field is induced in 
the neighbourhood of the body by magnetic field. The 
neutron stars and pulsars are examples of such bod- 
ies. The first model of electric field which is generated 
in the neighbourhood of a neutron star was developed 
by |Deutsch| (fl955|) . Some other models were suggested 



and studied by several authors ( Michel|1991 ) . All these 
models are based on suggestion that the neutron star 
is a conducting body. 

Conducting body rotating in its own magnetic field, 
generates a potential difference between the poles and 
equator. This potential difference leads to development 
of a co-rotating magnetosphere. By assuming that the 
magnetic field is dipolar, and unaffected by the trapped 
particles in the magnetosphere, and that the field dipole 
axis is parallel to the rotation axis, |Goldreich, P. fc Ju-| 



lian ( 1969 1 determined many of the properties of the 
magnetosphere. Electromagnetic field and magneto- 
sphere of the oblique rotator and conducting body have 
been studied by many authors. See for example [Cohen 



ney 



(1980) 



& Rosenblum (1972); Kaburaki (1981); Cohen & Kear 



Electromagnetic field of conducting body differs es- 
sentially from the pure dipole field. But there are also 
celestial bodies, which consist of non-conducting mat- 
ter and their axis of rotation is inclined with respect to 
the magnetic field axis. Magnetic field of such objects 
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in good approximation can be described as the field of 
an inclined rotating magnetic moment or "oblique ro- 
tator" (Babcock & Cowling 1953). Eelectromagnetic 



field in the near and intermediate zones of a magnetic 
dipole rotating in free space has been studied recently 



by Sarychev (2009). He has also estimated the energy 



which can be acquired by the particles during acceler- 
ation. 



In our previous paper (Epp & Masterova 2012 1 we where 



have calculated the field of rotating non-conducting 
body with the dipole magnetic field, and the dipole axis 
inclined to the rotation axis. An integral of motion was 
found and appropriate effective potential energy was 
studied. It was shown that each stationary point of the 
effective potential energy corresponds to particular so- 
lution of equations of motion for a charged particle. In 
this paper we develop investigation of the effective po- 
tential energy. We study here stability of the particle 
motion in the vicinity of the stationary points. Special 
attention is paid to investigation of regions accessible 
and inaccessible for a charged particle motion. These 
regions are described and depicted for particles with 
different values of the constant of motion. It is shown 
that there are closed allowed regions at certain values of 
the constant of motion. These regions are co-rotating 
with the electromagnetic field of the magnetized body. 



2 The field of precessing magnetic dipole 

Let us consider the field which is produced by precess- 
ing magnetic dipole. We define the law of motion of a 
dipole moment (J, in the Cartesian coordinate system as 
follows: 



[i = fi(shi a cos ujt, sin a sin u>t, cos a), 



(1) 



where /i > and oj > are the magnitude and angular 
velocity of the dipole, respectively, < a < ir is the 
angle between the vector fi and rotational axis. Con- 
sider the general formulas for the field of an arbitrary 
changing dipole ( Feynman||1964 ): 



(n x jii) (n x /i) 



(2) 



(n x (n x flj) 3n(n/^) - ii 3n(n/x) - /x 

" = 2 1 2 1 3 ' \ 6 > 

rc r^c r a 

where c is the speed of light, n = r/r is the unit vector, 
r is the radius- vector. Field is calculated at time t, and 
all the quantities on the right side of these equations 

T 

are calculated at the retarded time t = t . 

c 



We use a spherical coordinate system (r,9,tp), in 
which the components of electric field vector have the 
form: 



r 3 E r = 0, 

r 3 Eg — fip sin ql (p sin t — cos t) , 



(4) 
(5) 



r 3 E v = — ppcos#sina (pcosr + sinr) , (6) 



T = Wt'-if, p = — . (7) 
c 

The magnetic field vector has components: 

r 3 H r = 2p [sin a sin6>(cos r — p sinr) 

+ cos 9 cos a] , (8) 
r 3 Hg = — p [cos 9 sin a(cosr — p sinr — p 2 cost) 

— sin 9 cos a], (9) 
r 3 i?0 = —/j sin a(sinr + p cost — p 2 sinr) . (10) 

The same formulas for a — 7r/2 have been obtained by 
Sarychev|(2009[ ). It follows from Eqs Q - ^ that the 
components of electric field satisfy equation: 



Eg_ 
Eq8 



where Ei 



Eu 



E, 



= 1, 



p[i sm a 



(11) 



E 0ip — E a g cos 9. Hence, 



the vector E circumscribes an ellipse with semiaxises 
Eq0 and Eq v in the plane orthogonal to the radius vec- 
tor. This ellipse degenerates to a circle in the direction 
of axis of rotation (9 = 0, 7r), but in equatorial plane 

7T 

9 = — vector E oscillates in meridian plane (parallel to 
z-axis). 

It follows from Eqs ^ and ^ that the magnetic 
field has constant components (in contrast to the elec- 
tric field): 



2/x 



H rc = -j cos t) cos a, 



Hf)r = 



sm v cos a. 



These components relate to the field of a static mag- 
netic dipole if the inclination angle a is equal to zero. 
The electric field tends to zero in this case. 

One can see that time dependence appears in the 
formulas for electric and magnetic fields as composi- 
tion ujf — (p. It means that change of uit' is equivalent 
to change of ip. In other words, geometry of the electric 
and magnetic field looks like as the field rotates with an- 
gular velocity lo around z-axis. This conclusion relates 
also to geometry of the field lines of electric and mag- 
netic fields. At first glance it would seem that we have 
a paradox: the linear velocity of motion is getting more 
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then speed of light well away from z-axis. But motion 
of the held lines does not relate to transfer of matter or 
held energy. The above equations state only that the 
electromagnetic field at any point (r, 9, ip) of space is 
equal to value of the held at the point (r, 9,<p — Sip) at 
the moment t — Sip/uj. Evidently, in the far field zone 
only radiation field remains, which moves radially with 
the speed of light. 



3 Effective potential energy 

3.1 Dynamics of a charged particle 

Let us hnd the Lagrangian for a charged particle in the 



held of precessing magnetic dipole ( jLandau k, Lifshitz 
19751): 



1 - 



{Av). 



(12) 



where v is the particle velocity vector, A is the vector 
potential. It is easy to prove that the vector potential 



(n x /j.) (n x fi) 



(13) 



gives the field ^ and Substituting spherical com- 
ponents of vectors A and v 



j4 ,-• 
A 

Am 



= 0, 

fi sin a 



(sinr + pcosr), 



fi 



cos a sin ( 



— sin a cos 9 (cos r — p sin r) 
v = (f, 9r, iprshid) 



(14) 
(15) 

(16) 
(17) 



into Lagrangian (12), we obtain 



ip 2 r 2 sin# 2 



cr 



1 9 sin a(sin t + p cos r) + tp sin 2 9 cos a 



-<^> sin # cos # sin a(cosr — psinr)} . 



(18) 



This Lagrangian is a function of all spherical coordi- 
nates and time. The time dependence appears in this 
formula as combination cot' — ip. 

Substituting this function into Lagrange's equation 



d dL dL 
dt dq dq 



= 0, 



(19) 



we obtain the relativistic equations of motion of a 
charged particle in the field of precessing magnetic 
dipole (see Appendix [A|. 



Further we consider the charged particle to be a non- 
relativistic one. Let us introduce a new set of general- 
ized coordinates p, 9, ip with 



ip = (p — Lot. 



(20) 



Actually, this means that we use a co-rotating reference 
system. It is known that the rotating reference system 
is applicable only inside the light cylinder, i.e. rw < c 
or p < 1. Then the nonrelativistic Lagrangian function 
takes the form 



L 



mc 
2^2 

c 2 p 



p 2 + p 2 9 2 + p 2 (ijj + u) 2 sin 2 9 (21) 
sin a (ip +w)sm0cos0(cos£ + /?sin£) 

+ 0(sin£ — pcos^) + e ^ J - cos a(ip + uj) sin 2 9, 
i c 1 p 

where £ = ip+p. Now, the Lagrangian and Hamiltonian 
functions do not depend on time explicitly. In this case 
the value of the Hamiltonian is an integral of motion 



H = 



dL mc 2 , 9 

Qanr - L = tt^O 9 



' dq 
1 
2 

e/xw 2 sin 29 



2lo 2 



2 2 • 2 

mc p sin I 



efiuj 
c 2 p 



o 2 9 2 +p 2 sin 2 9^P 2 ) 



cos a sin 2 9 



2c 2 p 



sina(cos£ + psin£). 



(22) 



The first term in this expression K = — - (p + p~ 

p 2 sin 2 9tp 2 ) is always positive and it can be considered 
as the kinetic energy. The rest terms are usually re- 
ferred to as effective potential energy. It can be written 
as: 



N i sin 26* 



ef 




(cos£ + psin£) 



(23) 



N, = 



e^Liu sin a 



Nn = 



epuj cos a 



In this notation formula ( 22 ) can be represented as 



K = H-V e 



(24) 



The inequality % — V e f > imposes restrictions on 
possible area of the particle motion. 

3.2 Stationary points of the effective potential energy 

Let us investigate the effective potential energy. In the 
stationary points of the effective potential energy the 
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particle can be in a stable, unstable or indifferent equi- 
librium. The aim of this investigation is the stability 
of a charged particle at the stationary points. As the 
Lagrangian is written for a non-relativistic particle, we 
consider that cur <C c or p <C 1. This means that the 
particles moving around the axis of precession with an 
angular velocity of about oj are non-relativistic. In this 
approximation f w ip and the effective potential energy 
takes the form 



ef 



-p 2 sin 2 9 



Nj_ sin 26» 



COS Ip 



2N» sin 2 6 



(25) 



In order to obtain the stationary points of the function 
V e f we find solutions of the set of equations: 



dV f 



ef 



0. 



(26) 



where = p,6,tp. Hence, we have a system of three 
equations 



- 2p 3 sin 2 9 + 27V|| sin 2 6 - Nj_ sin 26 cos ip = 0, 
-p 3 sin 26* - 2iVj| sin 26 + 2N±_ cos 26 cos V = 0, 

Nj_ sin 20 sin ^> = 0. 



Equation (29) has next solutions: 



Ip = 0, 7T, 

(n e Z). 



(27) 
(28) 
(29) 



(30) 
(31) 



The sign before the square root is defined by the sign 
of the charge which is hidden in N, and by condition 



p > 0. Therefore, Eq. (35) takes the form 

<. N 



cos a + q 



sin 2 a 



(36) 



where q = e/|e| = ±1 is the sign of the charge. Inserting 



p into Eq. (32), we find: 

r — 3 cos a + q\f 9 — sin 2 i 

2 sin a L 



tan 6> 



(37) 



Thus, considering (30), (35), (37), for the positive 



particles in the case when sin 26 ^ we get the follow- 
ing series of coordinates of stationary points for V e f ■ 



Pi 



N 



4 

tan 6*i = — 



cos a + 
1 



sin 2 a 



2 sin a 



3 cos a + \/ 9 — sin 2 i 



(38) 



N 



r>2 



4 

tan 62 = 



cos a + V9 — 
1 



sin 2 a 



i'2 = 7T, 



2 sin a 



3 cos a + V» — sin 2 1 



(39) 



For the negatively charged particles, we have next two 
stationary points: 



4 



N r 



tan 8 3 



cos a 



1 

2 sin a 



V9~ 



sin 2 a 



^3=0, 



-3 cos a + v9 — sin 2 . 



(40) 



3.2.1 Solution for ij) — 0, 7r 



Using Eq. ( 30 1 we can eliminate the variable tp from 



the equations (27) and (28) by substitution cos^ = 



e, where e = +1 corresponds to ip = and e = — 1 
corresponds to ip = 7T. This results in a system of two 
equations with two unknowns: 

- p 3 + N\\ - eN ± cot 6 = (32) 
-p 3 cot 6> - 27V,, cot 6 + eiVj_ (cot 2 9 - 1) = (33) 



pi 



2V 



4 

tan 04 



cos a 



^9~ 



sin 2 a 



Ip4 



' r 3 cos a - V9 



2 sin a 



sin a 



(41) 



Hence, the solution of equations (27) - (29) for the 



case sin 29 ^ gives two stationary points for a positive 
charge and two points for a negative charge. 

'KTL 

3.2.2 Solution for 9 = — 



Expressing cot 9 from equation (|32j) and substituting it 

(34) 



into equation ( 33 ) , we obtain 

2p 6 - N n p 3 - N 2 = 0, 



where N 

3 N 
P = — 
H 4 



epuj 
mc 4 



This equation has a solution: 



cos a 



± v/sT 



sin 2 a 



(35) 



It follows from Eqs (27) - (29) that at the axis of ro- 
tation (9 = 0, 7r) all the first derivatives from effective 
potential energy are equal to zero only in the plane 
cos ip = and for any p. Which means that on the axis 

9 = 0, 7r there are not stationary points. 

7T 1 1 1 

As to the equatorial plane 9 = —, Eqs (27) - (29) 

have the next solution: 



p = N,f , cosip = 0, iVii > 0. 



(42) 
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This gives two solutions for ip, because ip £ (0, 2n). 

7T 

Thus, in case 9 — — there are two stationary points: 



P5 



P6 



2' ^ 



71" 

= 2' 

T' 



(43) 
(44) 



It follows from (42) that ecosa > 0, which means that 



the two above mentioned stationary points correspond 
to a positive charge if a < tt/2 and to a negative charge 
if a > tt/2. 

Particle located at a stationary point can be in a 
state of equilibrium. Let us verify whether a particle 



with initial coordinates (pi,9i,ipi) defined by Eqs (38) 



(44) and initial zero velocity with respect to the ro- 



tating reference frame will be in equilibrium position. 
Substituting these coordinates and p = 9 = ip = in 



equations of motion (All - (A3) and taking into ac- 
count that ip = cut + ip we obtain identical equalities. 
Hence, a particle being at rest at one of the stationary 
points in the co-rotating coordinate system will keep 
this position. This means that in laboratory reference 
frame the particle is moving along a circle with con- 
stant velocity Vi = cpt sin 9i. Thus, there are at least 
six particular solutions which describe circular motion 
of the particles in the field of inclined rotating dipole. 
Positions of the orbits defined by angle 9 and their ra- 
dius are different for the positive and negative parti- 
cle. Two of the trajectories are laying in the equatorial 
plane z = 0. 

3.3 Labelling the stationary points 



In order to define whether the stationary points ( 38 ) 



(44) represent maximum or minimum of potential 



energy we apply Sylvester's criterion: suppose that in 
some neighbourhood of a stationary point Mi{pi, 
the second derivatives of the function V e f(p,9,ip) are 
defined. The function has a local maximum if the 
quadratic form d 2 V e f(Mo) is negative definite and a 
local minimum if d 2 V e f(Mi) is positive definite. If the 
form d 2 f(Mi) is alternating then there is no extremum 
at the point Mi (Merkin 1996). Let us find all sec- 



ond partial derivatives and define their values at the 



stationary points (38) - (41) 
2 d 2 V 



an 

«22 

033 

0,12 

a 2 3 

«13 



mc 2 dp 2 
—q cos av9 

2 d 2 V 
mc 2 d9 2 

2 d 2 V 



cos a 



sin 2 a, 



6p 2 



N sin 2 a 



mc 2 dtp 2 3p 
—qVd — sin 2 a), 



- (cos a 



2 d 2 V 
mc 2 dpd9 

2 d 2 V 
mc 2 d9dip 

2 3 2 V 
mc 2 dpdip 



2N cos ip sin a 

^2 ~ ' 



= o, 

= 0. 



(45) 
(46) 

(47) 
(48) 

(49) 

(50) 



Using Eq. (37) we transform Eq. (45) into a n 



-6 sin 9 and hence, it takes only negative values. 
For the first fundamental form, we have: 

2N r „ 



Di = 



an a 12 

a 21 a 22 



—a cos a 



cos 3 a 



■ 2 

sin a 



+ (cos 2 a - 4)gV9~ 
The last equation can be transformed to 



(51) 



P 



sin 2 a, 



which makes obvious that D\ < both for negative and 
positive particles. The third order determinant 



an 0,12 ai3 

Ct21 022 &23 
O31 032 »33 



= -E>ia33 



(52) 



is definitely negative because 033 > 0. Thus, the form 
d 2 f(Mi) is alternating and there are no extremums of 
the function V e f at the four points with coordinates 

(38b - m . 



Calculation of the second derivatives in case of 9 
gives: 



7T 

2 



an 



-6, a 22 = QN.f , a 33 = 0, 



aw = 0, (I13 = 0, a 2 3 = 2^^^!! 3 sin ip. 

The quadratic form combined of these derivatives is in- 
definite in its sign. 

Summing up, we state that the six stationary points 



( 38 ) - ( 41 ) , ( 43 ) and ( 44 ) are points of equilibrium for a 



particle at rest, but the equilibrium is not stable. "Par- 
ticle at rest" means that coordinates p, 9, ip are fixed. 
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In the inertial reference frame the particle is moving in 
a circle with constant velocity. 

Stationary points that we considered above do not 
represent all solutions of the equations of motion for 
circular orbits. If we substitute into equations of motion 
(All - (A3) p — const we find one more solution for 



ecosa > 0: 
2' F 



2N, ip = nt + Lp , 



(53) 



where f2 = (l/2)w, and (po is an arbitrary constant. 
This solution is valid for a positive charge if a < ir/2 
and for a negative one if a > ir/2. 



4 Equipotential surfaces of the effective 
potential energy 



Let us consider a case of positively charged particle. 
(56) has not positive roots for p if C > and 
inequality ( 55 ) is satisfied in all space. There 



Eq 



is no 



inaccessible area for particle motion. 

In case C £ (— oo; —3) Eq. (56) has two roots at any 



9. The typical form of the equipotential surface in this 
case is given in Fi g. [I] The inner surface has form of a 
torus. Inequality (|55| is satisfied inside the torus - this 
is a region accessible for motion. The region between 
the torus and outer surface is inaccessible. Outside the 
external surface we have again accessible region. As the 
constant C increases from— oo to —3, the inner surface 
of the torus is increasing, while the radius of outer sur- 
face is decreasing. At C = —3 the surfaces adjoin to 
each other along a circle of radius p — 1 . This is shown 
in Fig. [2] 



Let us return to the expression ( 24 ) : K = U— V e f ■ The 
integral of motion H is defined by the initial coordinates 
and velocity of the particle. If these are specified, the 
particle can move only in space area where V e f < H. 
In this section we present plots of surfaces V e f = const. 
These surfaces restrict the permissible area for particle 
motion. 

As we found in the previous section, the stationary 
points arc at distance p of order | | 1 / 3 from the co- 
ordinate origin. Hence, we introduce a variable p = 
p\N\~ x / 3 in order to make p ~ 1, and a reduced poten- 
tial energy V = 2V ef /(mc 2 A 2 / 3 ) 




Fig. 1 Equipotential surfaces for a — 0, C — —3.5, q = 1. 
C and A are accessible regions, B - inaccessible region 



- _ 2 . 2 q sin a sin 29 

V = — p sin H cos?/' 



2q cos a sin"" 6 



(54) 



where q = e/|e|. The regions accessible to the particle 
motion are defined by inequality 



V < C. 



(55) 



The constant C is related to the integral of motion T-L 
as C = 2H/{mc 2 N 2 / 3 ). 

4.1 Aligned rotator, a = 

It is useful to start drawing plots from the simple case 
of a stationary dipole a = in order to show how the 
axial symmetric equipotential surfaces transform into 
antisymmetric ones. If a — 0, the equipotential surface 
for a given constant C is defined by equation (9 7^ 0, 7r): 



p 3 + T]p + 2q = 0, 
where r\ — C j sin 2 9. 



(56) 




Fig. 2 Equipotential surface for a = 0, C — —3, q = 1. A 
and C are accessible regions, B - inaccessible region. 



In interval C e (-3;0) Eq, 
roots if sin 2 9 < — C/3, one root if sin 2 9 



( 56 ) has two positive 
-C/3 and 



no roots in case sin 9 > —C/3 . The typical form of the 
equipotential surface in the case C € (—3; 0) is given in 
Fig. [3] The space inside the surface is inaccessible for 
motion, while the outside space is accessible. The upper 
and lower parts of the surface in Fig. [3] have point-like 
connection at the coordinate origin p = 0. In the limit 
C — > —0 all the space is accessible for motion. 
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Fig. 3 Equipotential surface for a = 0, C = 
is inaccessible region, A - accessible region 



-2, q = l. B 



For a negatively charged particle equation ( 56 1 has 



one root for each C. In the case C G (0, +oo) we have 
a surface in the form of a torus (Fig. Ej). In other cases 
(C G (0, — oo)), the equipotential surface has a form 
shown in Fig. [5j The radius of the surface increases 
with decreasing of C. 




Fig. 4 Equipotential surfaces at C = 3, q 
cessible region, A - accessible region. 



-1, B - inac- 



cquipotential surfaces are symmetric with respect to 
substitution 9 — > tt — 9, tj) — > ij) + tt, i.e. they are 
symmetric with respect to reflection in the plane z = 
and rotation around axis OZ through tt. 



Instead of ( 56 ) we have next equation 



-r]p + 2x = 0, X — <z(cosa — sin a cot cos ■0) . (57) 



In case x < this equation has one solution for p, hence, 
there is one equipotential surface which divides all space 
into accessible and inaccessible regions. If x > and 
C > the equation has no solutions, all space is acces- 
sible for the particle motion. If x > and C < 0, the 
number of solutions depends on the sig n o f expression 
6 = X — (— 77 / 3) 3 / 2 . Namely, equation (57 1 has no so- 
lutions for p if b > 0, one solution for 6 = and two 
solutions if b < 0. 

The sign of b is the same as the sign of expression 



q sin 2 #(cos a sin 9 — sin a cos 9 cos tp) 



-3/2 



Evidently, the sum of first two terms lies in the range 
[-1, 1]. Hence, if C < —3, b is definitely negative. In 
other cases sign of b depends on 9 and a. 





Fig. 5 Equipotential surface at C = — 2,q — — 1. 
accessible region, B - inaccessible region. 



A is 



4.2 General case, a ^ 

Let us consider the equipotential surfaces in the general 



case a =^ 0. One can see from equation (54) that the 



Fig. 6 Equipotential surfaces at C = —4, q ■ 
are accessible regions, B - inaccessible region. 



1. A and C 



Thus, if C € (— oo; —3), equation (57 1 has two roots 
at any 9. The typical form of the equipotential surface 
in this case is given in Fig. [6j But unlike the case a — 
(Fig. [I]), there are three equipotential surfaces. 

Equipotential surfaces for other values of C and q 
are shown in Figs [7]|9] All pictures |6][9] are plotted for 
a = 7r/3. 

4.3 Stationary points and equipotential surfaces 

Let us prove that the stationary points found in section 



3.2 are located in accessible regions. One can see that 
in the limit a — > all these points tend to points with 
coordinates 9 = tt/2, p = N 1 / 3 ,^ = 0, tt, tt/2, 3tt/2. A 
particle being at rest at these points in the rotating ref- 
erence frame, is moving along a circle in the laboratory 
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Fig. 7 Equipotential surface at C = — l,q = 
are accessible regions, B - inaccessible region. 



1. A and C 



Then, for the stationary points ( 38 1 , ( 39 ) we have C ps 
—2.314. The equipotential surface for this value of C 
is shown in Fig. 
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The stationary points are points 
of contiguity of regions A and C. Again, a particle at 
these two positions is in unstable state of motion. 

-0.816. 
We see that points 



For the solutions ( 40 ) , ( 41 ) we have C 
This case is depicted in Fig. |11| 



(40), (41) are points of unstable motion. 





Fig. 8 Equipotential surface at C = 3, q = — 1. 
accessible region, B and C - inaccessible region. 



A is 




Fig. 10 Equipotential surfaces at C — — 2.314, q = 1. 
and C are accessible regions, B - inaccessible region. 




Fig. 11 Equipotential surfaces at C = — 0.816, q = 
and C are accessible regions, B - inaccessible region. 



-1. A 



Fig. 9 Equipotential surfaces at C — — 3, q = — 1. A and 
C are accessible regions, B - inaccessible region. 



reference frame. Substituting the values of coordinates 
and velocity in Eq. ( [22] ) we find that C — —3. Hence, 
the particle moves along the border between two regions 
A and C of Fig(2] This orbit corresponds to unstable 
motion. 

Let us find the constant C for stationary points in 



case a^0. Substituting (36) and (37) in equation (54) 
we obtain: 



C = — 



^3 cos a + q\J 9 — sin 2 



2V3 



cos a + q 



sin 2 a 



1/3' 



(58) 



5 Conclusions 

We have calculated the components of the magnetic 
and electric fields of a rotating non-conducting body 
when the axis of the dipole magnetic field of the body 
is inclined to the rotation axis. 

In order to find whether the "radiation belts" are 
possible in such system, we have investigated the ef- 
fective potential energy in the field of inclined rotat- 
ing magnetic dipole. It is shown that there are closed 
equipotential surfaces which are co-rotating with the 
dipole field. The particles with definite initial energy 
are confined within such surfaces. The circular orbits 
of charged particles correspond to the points at which 
the two allowed for motion regions are contiguous to 
each other. It is shown that these orbits are unstable. 
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Equipotential surfaces are plotted for different val- 
ues of the constant of motion for positive and negative 
charge of a particle. The results can be used for descrip- 
tion of the radiation belts around some specific celestial 
bodies possessing inclined dipolar magnetic field. 

This research has been supported by the grant for 
LRSS, project No 224.2012.2 
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A Relativistic equations of motion 

Here we record relativistic equations of motion of a charged particle in the field of rotating magnetic dipole. 



Lagrangian (181 leads to the following equations of motion in a spherical coordinate system. 
rn 



(1-/32)3/2 y c 



r(c 2 



r 2 2 - r 2 tp 2 s i n 2 0) + — (f9 2 + r99 + ip sin 2 9 {rip + rip) 



1 



-^rip 2 9 sin 26>) } - 



9 2 + ip 2 sin 2 i 



ip sin 20(cos r — p sin r — p 2 cos r) } 



c 

e/j, sin a 
2cr 2 



1 26 '(sin t + p cos t — p 2 sin r) 



e/i cos a . . 2 

— -= — yjsur0 = 0, 



(Al) 



2(! _ £2)3/2 



r 2 9(c 2 - r 



2 r 2 2 sin 2 9) + rr9{2c 2 ~ 2f 2 - r 2 9 2 - r 2 p 2 sin 2 i 



o . .. . n . 1 9 .OA . „„, 1 mip 2 r 2 cos0sin0 
+r 2 9{rr + r 2 ipp sin 2 + -r 2 <i> 2 sin 20) ^ y 

2 J 

e/isina r . o/iw «. . . i such cos a sin 20 

7. \ p(uj — p — 2p sm 0) (cos t — p sin r ) — p sin r f = 0, 

r z LU L J cr 



(A2) 



(1 - /3 2 ) 3 / 2 



sin 2 8(2rf(p + r 2 p) + r 2 <^0 sin 20 



r sin 2 



r^(2f 2 + r 2 2 + rV sin 2 0) 



+r(,5(f 2 + r 2 9 2 ) - rip{rr + r 2 99 + -r 2 p 2 9 sin 20) 



- sin29(r 2 + r 2 9 2 + r 2 tp 2 sin 2 ■ 



pte sin a J 2r0 sin 



- (cos t — p sin r) 



f sin 20 



cos r — p sin r — p 2 cos r) 



psin20 / . | p,ecosa f ^ . . 



(sinr + pcosr) > + 



?r sin 20 - r sin 2 0) = 0. 



(A3) 



B Stability of the orbit in the equatorial plane 



It is easy to check that equations of motion ( Al ) - (A3 1 have partial solution 



r = R, 6 = -, ip = ujt + ip a 



(Bl) 



with restiction on the initial phase p>q: 

cos(po + Po) + Po sin(p + ip ) = 0, p Q = Rw/c. 

We assume now that the particle is a non relativistic one. Then p$ I & n d we have for the initial phase cos po = 0. 



This solution corresponds to motion of a particle being at rest at the stationary points ( 43 ) and ( 44 1 . Let us verify 
stability of this solution. We develop equations of motion as series in increments of coordinates in a neighbourhood 
of the solution. Let us introduce the increments of coordinates as follows 



= — + 69, p = cut + po + Sp, r = R + Sr. 



(B2) 



We substitute these coordinates in equations (Al ) - (A3) and expand the equations up to first order of coordinates 



increments. Taking into account that for a non reativistic particle p <C 1 we obtain a system of differential equations 
that describe the motion of a particle in a vicinity of a circle: 



Sr — 3Sruj — uiRSip + eujR tg a& 
Sfutga , a 2 



0. 



R 



iuj 2 89~euj z tgaSp = 0, 
Stuj 

Sp + — — - ecu 2 tg aS9 = 0, 
R 



(B3) 
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where e = sin (po = ±1. We investigate the trivial solution of this system for stability. By Lyapunov's theorem the 
trivial solution is stable when all roots of characteristic first-approximation equation of the system have negative 
real parts ( Merkin|[l996 ) . We write Eq. (B3| in the form: 



— 3xi — x' 3 + Tx' 2 = 0, 
Tx\ + 3x2 — Tx 3 = 0, 
x'i + x\ - Tx 2 = 0. 



(B4) 



where x\ 



R 



,x 2 



,X3 = Sip,T — etga. Prime denotes £ derivative , where £ = lot. Reduce the system (B4) 



to a normal form: 



x\ - di = d[- 3xi - d 3 + Td 2 = 
d' 2 — Td\ + 3x 2 — TX3 = 
d' 3 + di- Tx 2 = 



x' 2 — d 2 = 
x' 3 -d 3 = 



(B5) 



A set of solutions of the system has the form d\ — aie w ^,d 2 — a 2 e lv( - ....x 3 — a^e w ^. This gives the characteristic 
equation: 



k 3 - fc 2 (l + T 2 ) - 3fc(l + T 2 ) - 3T 2 = 0, 



(B6) 



where k — v 2 . 

Important criteria that give necessary and sufficient conditions for all the roots of the characteristic polynomial 



(with real coefficients) to lie in the left half of the complex plane are known as Routh-Hurwitz criteria ( Gantmacher 
[I959] ). 

Given a polynomial 



P(A) = A™ + aiA"- 1 + ... + a„_!A + a n , 



(B7) 



where the coefficients are real constants, i = 1, 2, 3..., n. Define the n Hurwitz matrices using the coefficients 
of the characteristic polynomial: 



H, = (ai), H 2 



ai 1 

«3 a 2 



, H 3 




(B8) 



All of the roots of the polynomial P(A) are negative or have negative real part if the determinants of all Hurwitz 



matrices are positive (DetiJj > 0, j = 1, 2, . . . , n). 

For the polynomials of degree n=3, the Routh-Hurwitz criteria simplify to 

ai > 0, a\a 2 — a 3 > 0, a 3 (a\a 2 — a 3 ) > 0. 



(B9) 



From (|B6| wc get: a a = -(1 + T 2 ), a 2 = -3(1 + T 2 ), a 3 = -3T 2 . It is easy to see that conditions (|B9|) of 
Routh-Hurwitz criteria are not satisfied. So, the trivial solution of system (B3) is not stable. 
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